The piston dispersive shock wave problem 
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The one-dimensional piston shock problem is a classical result of shock wave theory. In this 
work, the analogous dispersive shock wave (DSW) problem for a dispersive fluid described by the 
nonlinear Schrodinger equation is analyzed. Asymptotic solutions are calculated using Whitham 
averaging theory for a "piston" (step potential) moving with uniform speed into a dispersive fluid 
at rest. These asymptotic results agree quantitatively with numerical simulations. It is shown that 
the behavior of these solutions is quite different from their classical counterparts. In particular, 
the shock structure depends on the speed of the piston. These results have direct application to 
Bose- Einstein condensates and the propagation of light through a nonlinear, defocusing medium. 
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The study of dispersive shock waves (DSWs) has 
gained interest with the recent experimental reaUzation 
of DSWs in a Bose-Einstein condensate (BEC) and 
the propagation of Ught through a nonhnear, defocus- 
ing medium Comparisons between classical, viscous 
shock waves (VSWs) and DSWs have been discussed in 
the context of single shocks [l[ and the interaction of 
two shocks Q yielding some appealing similarities but 
also important differences. Motivated by the classical 
VSW piston problem, here we consider the generation of 
a DSW by a piston moving into a dispersive fluid at rest. 

The theoretical study of DSWs involves averaging a 
periodic wave over its period and allowing for slow vari- 
ation of the wave's parameters. This method, known 
as Whitham averaging has been successfully applied 
to many DSW problems including step initial data for 
the nonlinear Schrodinger (NLS) equation 0, 0] , Bose- 
Einstein condensates EI.Bl , fiber optics , the generation 
of ultrashort lasers [To], and DSW interactions Qj- We 
also note that a dispersive piston shock problem arises 
as an asymptotic reduction of 2D, supersonic flow of a 
dispersive fluid around an obstacle 

The piston shock problem is one of the canonical prob- 
lems in the theory of VSWs. A uniform gas is held at 
rest in a long, cylindrical chamber with a piston at one 
end. When the piston is impulsively moved into the gas 
with constant speed, a region of higher density builds 
up between the piston and a shock front which propa- 
gates ahead of it. An elegant asymptotic solution to this 
problem is well known and relates the shock speed to the 
speed of the piston and the initial density of the gas (see 
e.g. and the discussion at the end of this work). 

In this work, we consider the problem of a "piston" 
moving with constant speed into a steady, dispersive 
fluid: e. g. a Bose-Einstein condensate or light propagat- 
ing through a nonlinear, defocusing medium. The piston 
in this case is a step potential that moves with uniform 
velocity. This potential could be realized in a BEC with 



a repulsive dipole beam and in nonlinear optics with a 
local change in the index of refraction. One expects, in 
analogy with the classical, viscous case, the generation 
of a dispersive shock wave. As we will show, this is in- 
deed the case. There are two types of asymptotic behav- 
ior depending on the piston velocity. For smaller piston 
velocities, a region of larger density or intensity builds 
up between the piston and a DSW. However, for large 
enough piston velocities, a locally periodic wave is gen- 
erated between the piston and the DSW which has no 
VSW correlate. The asymptotic results are verified by 
numerical simulations demonstrating quantitative agree- 
ment. 

We consider the ID NLS equation with a potential 
(also known as the Gross-Pitaevskii (CP) equation) 



-^.^ + Voix,t)^ + 0<e«l. (1) 



This equation models the mean field of a quasi-lD BEC 
[ist and the slowly varying envelope of the electro- 
magnetic field propagating through a Kerr medium [l3 | 
(where time t is replaced by the propagation distance). 
The small parameter e is inversely proportional to the 
number of atoms in the BEC [H or, after rescaling, in- 
versely proportional to the maximum initial intensity of 
the electromagnetic field. For all calculations in this 
work, we assume e = 0.015, a typical experimental value 
for BEC [ij . The piston problem is modeled with a tem- 
porally and spatially varying step potential given by 



Vb(x, t) = VmaxHiVpt ~ x), H{y) 



?/ < 
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with strength Vmax and constant velocity Vp. The initial 
conditions are 



^{x,Q) 



fpp^ as X ^ oo, ^(x, 0) ^ as a; — oo. 



Because the strength of the piston is large, Vmax ^ P_r, 
the density /intensity rapidly decay to zero near the ori- 
gin. We assume that the wavefunction is in the ground 



2 



state of the step potential VmaxH{—x) when t < 0. For 
all calculations in this work, pji = 0.133. 

It is useful to view eq. ([1]) in its hydrodynamic form 
by making the transformation \1/ = y^pe^" ■'o '^^ and 

inserting this expression into the first two local conserva- 
tion equations for the GP equation 



Pt + {pu)x = 



{pu)t + {pu^ + ip^)^ = — (p(logp)^^)^ - pVq^ 



where p is the dispersive fluid "density" and u is the 
dispersive fluid "velocity" . These equations are similar 
to the Navier-Stokes (NS) equations of fluid dynamics 
except that the viscous term of NS has been replaced by 
the dispersive term with coefficient 

Because the dispersive term in eq. ^ is small, one ex- 
pects the generation of small wavelength 0(e) oscillations 
near a steep gradient in the fiuid variables. Witham's 
method is to average an exact periodic solution over fast 
oscillations and assume that the wave's parameters (am- 
plitude, frequency, wavelength, etc.) vary slowly 0. 

We convert the piston DSW problem into a moving 
boundary value problem where appropriate boundary 
conditions are imposed at the piston front. First we solve 
the piston DSW problem assuming sufficiently small, pos- 
itive piston velocities Vp. "Small" will be defined below. 

We assume the piston strength is large for x < Vpt, 
so there is negligible density there. Assuming there is a 
jump from zero density to the nonzero value pL with a 
fluid velocity wc integrate the first conservation law 
in eqs. ([2]) across the jump to flnd 

- VpPL + Plul = 0. (3) 
This gives the first boundary condition at the piston 



u{vpt, t) = ul = Vp, 



(4) 



the fluid velocity at the piston equals the piston velocity. 
We require a boundary condition for the density. 

The theory of DSWs involves a system of quasi-linear, 
first order, hyperbolic equations known as the Whitham 
modulation equations. The Whitham equations describe 
the slow evolution of a periodic wave's parameters and 
must be solved in order to find the asymptotic DSW so- 
lution. The simplest, non-trivial solutions to these equa- 
tions are known as simple waves where only one depen- 
dent variable is varying in space and time, and the rest 
are constant. In analogy with gas dynamics, we assume 
a simple wave solution, but in this case to the Whitham 
equations. This determines a density p^ at the piston. In 
order to connect to the uniform state ahead of the piston 
PR < PL, wc must have a single DSW for Vp sufficiently 
small (vp < 2y/pfi). As wc will show below, a "vacuum 
state" is created when Vp > 2y/pii, and we find a uniform 



traveling wave (TW) with speed Vp, instead of the con- 
stant density pi, adjacent to the DSW. Later we verify 
with numerical simulations that these assumptions are 
reasonable. Now we derive the asymptotic piston DSW. 

At the time i = 0+, we assume that there is a discon- 
tinuity in the fiuid variables due to the impulsive motion 
of the piston at t = 



(2) p{x,0+) 



PL X = _ \ UL =Vp X ^0 



PR X >Q 



, u{x,0+) 



ur 



x>Q ' 
(5) 

This discontinuity is regularized by a slowly modulated, 
traveling wave, periodic solution to eq. ([2]) with V{x, t) = 

i 



p{x, t, 6*) = A3 - [A3 - Ai]dn^(6'; to), to 



A2 — Ai 



/ , n\ M \/AlA2A3 



p{x,t,9) 



A3 — Ai 
±1, < Ai < A2 < A3 



de 



de 



dx v/^^/^' dt 



-y - Ai/£, 



(6) 



where the parameters A;, i = 1, 2, 3 and V satisfy 

Ai =11(^1 - ^2 - T-s + Tif, A2 = j^i-ri + r2 - r3 + r^f, 
A3 =l^(-^i - 7-2 + rs + r4)^ V = i(ri + ra + r3 + n), 

(7) 

and Ti = ri(x,t) evolve according to the Whitham equa- 
tions 



0, 



1,2,3,4. 



The velocities Vi are expressed in terms of complete el- 
liptic integrals of the first and second kind Q . 

In order to find a simple wave solution to the Whitham 
equations, we require that only one of the parameters 
ri spatially varies and that the initial data for all the 
parameters properly characterizes the initial data in 
eq. ^ with the spatial average of eq. ^ . We use the 
method of initial data regularization [j, [^, [lo| to find 



-"^y/PR., r2 = 2^pfl, r4 = 2vp + 2^/pfl, 



7-3 (x,0^ 



2Vp^ 







2vp + 2ypi a; > ' 



(8) 



p{vpt,t) = PL = ii:Vp + y/p^f. 



The last equation, the boundary condition for the density 
at the piston, comes from the simple wave assumption. 
Equations ^ give rise to a self similar solution for 
satisfying the implicit relation 



V3iri,r2,r3{x,t),r4) x/t. 



(9) 



Using a nonlinear root finder, eq. Q is solved numerically 
for each x and t. The values for r^, i = 1,2,3,4 are 
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FIG. 1: LEFT: Asymptotic piston DSW solution at time t = 3 
for Vp — y/pR = 0.365. The vertical lines mark the left and 
right edges of the DSW moving with speeds v~ = 0.548 and 
v+ = 1.278 respectively. RIGHT: Asymptotic piston DSW 
solution at time t — 1.7 when Vp — 2.5y^pij = 0.912. A locally 
periodic region connects the piston to the trailing edge of the 
DSW. The density minima in this region are approximately 
zero and the velocity is, theoretically, undefined (infinite) at 
these points. The density maxima are 4pR = 0.532. DSW 
speeds are = 1.245, vf = 2.452. 



inserted into eqs. ([7]) and ^ to determine the asymptotic 
DSW solution. A pure DSW propagates ahead of the 
piston with traiUng and leading edge speeds respectively 

m 



1 



/PR, 



+ 



PR + PR 



(10) 



Figure [TJ left depicts the asymptotic piston DSW solu- 
tion for a small piston velocity. The minimum values of 
the density and velocity occur at the trailing edge of the 
DSW and are jH, 



PR 



/PR 



P^)- (11) 



The maximum values occur between the piston and the 



DSW: p„ 



PL = ivp/2 + ^/prY, Ur. 



UL 



It is possible for the piston velocity to be greater than 
the trailing DSW velocity calculated using cq. pO|) 



if 



> 2 



When Vp = p can vanish (there is a so-called vac- 

uum point) and a modification of the solution is required. 
To find a simple wave solution for large piston velocities, 
we must derive new conditions for the parameters r,; . We 
modify the DSW solution by introducing a locally pe- 
riodic TW between the piston and the trailing edge of 
the DSW. When Vp ~ 2^fp^, the DSW forms a vacuum 
point [H, 0] at the piston. The vacuum condition, p = 0, 
is satisfied in cq. ^ when 

Ai = =^ ri-r2-r3+r4 = ^ r^-r^ ^ iy^. (12) 



FIG. 2: Numerical solutions to eq. ((ij for the piston problem. 
LEFT: The density (upper) and velocity (lower) for the same 
parameters as in Fig.[T]LEFT. Numerically calculated trailing 
edge speed is 0.557, approximately the theoretical value 0.548. 
RIGHT: The density (upper) and velocity (lower) for the same 
parameters as in Fig. [T] RIGHT. Numerically calculated TW 
velocity of locally periodic wave is 0.912, approximately the 
piston speed 0.913. 



Note that the fluid velocity is undefined at a vacuum 
point, even though the vacuum points have a well-defined 
propagation speed through the fluid. We assume that 
this condition holds for Vp > l^fJT^ as well. One more 
condition is required to completely determine r/^ and 
T-i(x = 0,t = 0+); ri, r2, and rj,{x > 0,0+) are deter- 
mined by the initial data ([5]). Because there is a locally 
periodic TW between the piston and DSW, we assume 
that the velocity of the TW equals the piston velocity. 
This is the TW velocity condition 



V 



ri+r2+r3+ri = 4wp => r^+r^ = Avp. (13) 



Given the initial data in eq. (O and the two conditions 
([T^ and p^ . only and a in the initial data of cq. ^ 
are altered 

r3{x,0+) =2vpT2,/^„ (7(2;, 0+) = Tl, x = 0. 

Note that a = —1 in the locally periodic region. Figure 
[U right depicts the asymptotic DSW solution for Vp > 

Several properties of this DSW solution are worth not- 
ing. The density between the piston and the DSW oscil- 
lates between the values 



Pr, 



and p„ 



4=PR, 



(14) 



independent of the piston velocity Vp and the TW in 
this region propagates with the velocity V = Vp. Note 
that since the vacuum condition (|12p holds everywhere 
inside the TW trailing the DSW, the velocity in this re- 



gion, from cq. 



V 



everywhere (except 



at the vacuum points where the velocity is undefined). 
The wavelength of the TW is I = ^K{4pii/v^), where 
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K{m) and E{m) are the complete elliptic integrals of the 
first and second kinds respectively. The DSW propagates 
with trailing edge speed (also the propagation speed of 
the rightmost vacuum point where tr changes sign) 



'"s = Vp + ivp + 3y/p]i) 



Vj,E{Apr/vI) 



{vp - 2^)K{ApR/vl) 



and leading edge speed v'^ , the same as that given in eq. 
(jlOp . The number of vacuum points increases linearly 
with time: 



I 



-t 



2eK{ApR/vl) 



(15) 



We perform direct numerical simulations of cq. ([T]) to 
verify the assumptions we have made such as the bound- 
ary conditions ([4|), ([8]), the vacuum and TW velocity con- 
ditions (fT^ . p^ . and the trailing edge DSW speed 
of eq. pO]) . All of our assumptions are in excellent agree- 
ment with numerical simulation as shown in Fig. [31 Nu- 
merical simulations of eq. ([T]) were performed with the 
pseudo-spectral, Fourier method [l^ with a grid spacing 
Aa; = 0.004, time step = 0.0006, parameter Vmax = 5, 
and a slightly smoothed step potential Vq. The initial 
data is relaxed in the presence of the potential Vq and 
spatially localized by including a smoothed step poten- 
tial with strength Vmax near the right boundary of the 
spatial domain. 

Numerically calculated piston DSWs for both moder- 
ate Vp = yfpR and large Vp = 2.byfpii piston velocities 
are shown in Fig. [51 For the slower piston velocity in 
Fig. [21 left, the solution is similar to the asymptotic re- 
sult in Fig. [H left. The only difference is the variation 
in the density and velocity through the piston. In the 
model problem considered here, we assume that the den- 
sity goes to zero immediately behind the piston. The 
piston DSW corresponding to the large piston velocity 
in Fig. [21 right is very similar to the asymptotic result in 
Fig. [H right. The vacuum condition in eq. (fT^ predicts 
u = Vp everywhere in the trailing wave region except at 
vacuum points where it is undefined. This is reflected 
in the numerical calculation as very large spikes in the 
velocity when the density approaches zero. 

The analogous piston viscous shock wave problem in 
shallow water is discussed in, e.g. the Id equations 
are equivalent to eqs. ([2]) when e = 0, F = and a 
dissipative regularization is used whenever a shock forms. 
The asymptotic solution is found by assuming a simple 
wave and incorporating the boundary condition In 
this case, one finds that the shock speed [vg) is always 
larger than the piston speed [vp), i.e. one finds Vg — Vp^ 
VpPR/{pL - Pr) > 0. 

The analysis in this work shows that techniques from 
VSW theory, simple wave solutions and suitable jump 
conditions, are useful in the study of DSWs. Neverthe- 
less, DSWs can lead to very different phenomena. 
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FIG. 3: Comparisons of the analytical theory (solid curves) 
and numerical simulations (dots) for different DSW parame- 
ters as functions of the piston velocity Vp. LEFT: pmax (top) 
and pmin (bottom). Maximum absolute error in pmax from 
eqs. (jSJ and (|14|l is 0.0088. Maximum absolute error in pmin 
from eq. (HTJ is 0.010. RIGHT: Top depicts the speed of the 
trailing edge of the piston DSW for Vp < 2ypij — 0.73 and 
TW velocity V of locally periodic wave for Vp > 2y/])R. Bot- 
tom shows the validity of the boundary condition Q. Maxi- 
mum absolute error in vl^ /V from eqs. (|10|l and (|13|l is 0.018. 
Maximum absolute error in ul from eq. Q is 0.013. 
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